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Abstract. This is the first of two articles deahng with the equa- 
tion {-AYv = f{v) in R", with s G (0,1), where (-A)" stands 
for the fractional Laplacian — the infinitesimal generator of a Levy 
process. This equation can be realized as a local linear degener- 
ate elliptic equation in M""*"^ together with a nonlinear Neumann 
boundary condition on 9R"+^ = R". 

In this first article, we establish necessary conditions on the 
nonlinearity / to admit certain type of solutions, with special in- 
terest in bounded increasing solutions in all of K. These necessary 
conditions (which will be proven in a follow-up paper to be also 
sufficient for the existence of a bounded increasing solution) are 
derived from an equality and an estimate involving a Hamiltonian 
— in the spirit of a result of Modica for the Laplacian. In addi- 
tion, we study regularity issues, as well as maximum and Harnack 
principles associated to the equation. 



1. Introduction 
This paper is devoted to the study of the nonhnear problem 

{-Ayv = f{v) inW, (1.1) 

where s G (0, 1) and 

(-A)^.(x) = c„,, p.v. / (1.2) 

Here P.V. stands for the Cauchy principal value, and C„ ^ is a positive 
constant depending only on n and s — whose value is given in Re- 
mark 13.111 below. The above integral is well defined if, for instance, v 
is bounded (which ensures the integrability at infinity) and v is Cj^^(M"') 
(which ensures the integrability at a; = x in the principal value sense). 
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As explained in more detail in section 3 below, up to an explicit 
multiplicative constant (given in Remark I3.1ip in front of the nonlin- 
earity /, problem (11 .ip can be realized in a local manner through the 
nonlinear boundary value problem 



div {y"- Vu) = in 

du ' ^1 (1.3) 

1 + a)— = f(u) on = M'^, 

where n > 1, = {{x,y) e x M : ?/ > 0} is the halfspace, 

= {?/ = 0}, M = ^(x, ?/) is real valued, and 

lim y°'dyU 



is the conormal exterior derivative of u. Points in M" are denoted by 
X = (xi, . . . , Xn)- The parameter a belongs to (—1, 1) and is related to 
the power s of the fractional Laplacian (—A)* by the formula 

a = 1 - 2s G (-1,1) 

— a relation that we assume throughout the paper. 

Indeed, Caffarelli and Silvestre [6] have proved the following formula 
relating the fractional Laplacian (— A)'^ to the Dirichlet to Neumann 
operator for (II ■3p : 

Hi I 

(-A)^{m(-, 0)} = d,— in = dWl+\ (1.4) 

where dg is a positive constant depending only on s (see Remark 13.111 
below for its value). 

The aim of this paper — and of the forthcoming one [1] — is to study 
two types of bounded solutions of (II. ip : 

(a) Solutions v = v{x) of (11.11) which are monotone increasing, say 
from —1 to 1, in one of the x- variables. These solutions are 
named layer solutions and constitute our main interest. 

(b) Radial solutions v = v{\x\) of (11.11) tending, say, to as |a;| — > 
±oo. 

In the second part [4| of this work, we will be concerned with the ex- 
istence, uniqueness, symmetry and variational properties, as well as the 
asymptotic behavior of layer solutions. These questions will be related, 
of course, to structural assumptions made on the nonlinearity /. 

In this first article, we establish necessary conditions on the nonlin- 
earity / to admit a layer solution in M (i.e., in dimension one), and 
also to admit radial solutions in M" with limit at infinity. In the case 
of layer solutions, our necessary conditions (12. 7p and (12.80 below will 
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be proven in [4] to be also sufficient for the existence of a layer solu- 
tion. Our necessary conditions are derived from a new equality and 
a new estimate involving the Hamiltonian for fll.Sp — in the spirit of 
a celebrated inequality of Modica [T6] for the Laplacian. In subsec- 
tion 1.1 below we explain the Hamiltonian structure of problem (11. 3p . 
Let us also recall that Modica proved that every bounded solution of 
Av - G'{v) = in the whole W satisfies (l/2)|Vt;p < G{v) pointwise 
in all M", assuming only that G > in M. Here we prove an analogue 
of such an estimate in dimension one in the case of nonlocal operators 
given by the fractional Laplacians. This is done via a careful study of 
the local boundary value problem ( II. 3p . 

In addition, in this ffist paper we also study regularity issues, as 
well as maximum, Liouville, and Harnack principles associated to the 
fractional Laplacian in M". These tools will be needed in this paper 
and in its follow-up. 

Our work extends to all fractions s G (0, 1) the results of one of 
the authors and J. Sola-Morales |5j for the case s = 1/2. [5J studies 
bounded (specially layer) solutions of 



which corresponds to the case a = in (II. 3p . It is well known that the 
Dirichlet to Neumann operator associated to the previous problem is 
precisely (— A)^/^. Therefore, layer solutions of the previous equation 
are actually heteroclinic connections (between —1 and 1) of 



where v is the trace of u on the boundary {y = 0} . The goal of our 
papers is to generalize this study to any fractional power s G (0, 1) 
of the Laplacian. We will make a great use of the tools developed 
in [5]. However, some new difficulties arise due to the degeneracy of 
the operator in (II. Sp . 

The first equation in (II. 3p is a linear degenerate elliptic equation with 
weight y°'. Since a G (—1, 1), the weight y"" belongs to the Muckenhoupt 
class see [17]. More precisely, a nonnegative function w defined in 
is said to be A2 if, for some constant C, 



for all balls 5 in M^. It is easy to verify that |?/|" G y42(M"+^) for 
a G (—1, 1). As we explain in section 3, this fact allows to use the 




fiy) inM^ 



n 
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regularity results of Fabes, Jerison, Kenig, and Serapioni [101 E] con- 
cerning divergence form equations with A2 weights. Another important 
property of the weight y"" is that it depends only on the extension vari- 
able y and not on the tangential variables x. Thus, (11. 3p is invariant 
under translations in x — as it is equation (11. ip . In [4J, this will allow 
us, for instance, to use for (II. 3p the sliding method of Berestycki and 
Nirenberg. 

The study of elliptic equations involving fractional powers of the 
Laplacian appears to be important in many physical situations in which 
one has to consider long-range or anomalous diffusions. From a proba- 
bilistic point of view, the fractional Laplacian is the infinitesimal gen- 
erator of a Levy process — see the book of Bertoin [3], for instance. 

A lot of interest is currently devoted to the study of nonlinear equa- 
tions involving fractional powers of the Laplacian. This type of op- 
erator seems to have a prevalent role in physical situations such as 
combustion (see |7]), dislocations in mechanical systems (see [13]) or in 
crystals (see [121 ED])- In addition, these operators arise in modelling 
diffusion and transport in a highly heterogeneous medium; they may 
be used in the description of the diffusion of a liquid in some heteroge- 
neous media, or as an effective diffusion in a limiting advection-diffusion 
equation with a random velocity field; see for instance [T5] . 

1.1. Formal Hamiltonian structure. As in [5], the quantity appear- 
ing in our main results — see for instance (12. 9 p below — arises naturally 
when one looks at problem (II. 3p for = 1 as a formal Hamiltonian 
system in infinite dimensions. Here the time variable is r = x, the 
position q is the function u{x, ■) = u{t, ■) in the halfline {y > 0}, and 
the momentum is p = q' = Ux{t, ■). 

From the action — that is, in PDF terminology the energy func- 
tional (12. ip below, which we already have at hand — we see that the 
Lagrangian is L{q,p) = (l/2)||p|||^ + W{q), with 



and \\w\\l^ = /q y''\w{y)\'^dy. 

The Legendre transform of L with respect to p gives the Hamiltonian 



Hiq,p) = / -{uUx,t)-ulix,t)}dt--—G{u{0,x)) 



m<l) = l\\dyq\\la + Y^G{qm 



1 




l\\p\\la-W{q). 
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One can easily check that its associated Hamiltonian system 




is formally problem (11 .Sp . 

Thus, our equation admits a Hamiltonian structure in an infinite 
dimensional phase space. However, in this paper we do not address the 
question of setting it as a true well posed semigroup. Note that a lot of 
challenging issues usually arise with infinite-dimensional Hamiltonian 
systems — see for instance [8]. 

1.2. Outline of the article. The paper is organized as follows. Sec- 
tion 2 contains the statements of our main results. In section 3 we 
explain the relation between problems (11. ip and (II. 3p . and we present 
the Poisson kernel and a regularity result for (II. Sp . Section 4 contains 
results on the operator La appearing in (II. 3p : we establish Schauder 
estimates, a Harnack inequality, a Hopf principle, maximum principles, 
and a Liouville theorem. Section 5 is concerned with the proof of the 
Hamiltonian equality and estimates. In section 6 we prove our results 
on layers as the fraction s tends to 1. Finally, in section 7 we collect 
the proofs of our main results. Theorems 12. 2| 12.31 12.41 and 12. 5[ using 
the results established in previous sections. 

2. Main Results: Hamiltonian identity and necessary 



Throughout the paper we assume that s G (0, 1) and that the non- 
linearity satisfies 



— which is defined up to an additive constant. We recall that the 
problem under study is (II. 3p . i.e.. 



conditions on the nonlinearity for existence 



/gC^'^(M) for some 7 > max(0, 1 - 2s). 
We will denote G the associated potential, i.e.. 



G' 



f 



div {y" Vu) = in R';+ 

(9?; 

(1 + a)— = /(n) onM: 



with a = 1 — 2s. In the sequel we will denote 

LaW = div {y" Vw). 
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We use the notation 

= {{x,y) eR''^' : y > 0,\{x,y)\ < R}, 

= {(x,0) G dRl+^ : \x\ < R}, and 

r+ = {{x,y)eR''+':y>0,\{x,y)\ = R}. 

It is easy to see that (11. 3p has a variational structure, corresponding 
to the energy functional 

E^+{u)= I -y''\Vu\'^dxdy+ [ -^—G{u)dx. (2.1) 
Jb+ 2 Jpo^ 1 + a 

We are concerned with the following types of solutions. The first 
class (layer solutions) consists of solutions which are increasing and 
have limits at infinity in one Euclidean variable in the space M" of x- 
variables. In the following definition, and for future convenience, after 
a rotation we may assume that such variable is the a;i-variable. 

Definition 2.1. We say that f is a layer solution of (II. ip if f is a 
solution of (II. ip satisfying 

v^, > in R" and 

lim v{x) = ±1 for every {x2, ■ ■ .,Xn) G ]R""\ (2.2) 

We say that m is a layer solution of (ll.3p if it is a solution of (II. 3p . 

u^, > on (9R!^_+\ and (2.3) 
lim u{x, 0) = ±1 for every {x2, . . . , x„) G M""^ (2.4) 

xi— >±oo 

It is important to emphasize that, for n > 2, the limits in (12.20 
and (12. 4p are taken for (x2, . . . ,Xn) fixed, and are not assumed to be 
uniform in (x2, . . . , x„) G M"~^. 

We will also study solutions v of (ll.ip which are radially symmetric 
(not necessarily decreasing) and such that 

lim v{\x\) = 0. (2.5) 

\x\—^ca 

We can now state our main results. The next theorem provides a 
necessary condition — (12.70 and (12.80 — on the nonlinearity / to admit 
a layer solution in M. In our subsequent paper [4|, this necessary con- 
dition will be proven to be also sufficient for the existence of a layer 
solution. It is interesting to point out that conditions (12. 7p and (12. 8p 
are independent of the fraction s G (0, 1), and that they are also the 
necessary and sufficient conditions for the existence of a layer solution 
to the local equation —v" = f{v) in all of M. 
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The theorem also states that famihes of layer solutions indexed by 
s G (0, 1) converge as s goes to 1 to a layer solution of the equation 
-v" = f {v) in R. 

Theorem 2.2. (i) Let s G (0, 1) and f any C^''^{W) function, for some 
7 > max(0, 1 — 2s). Assume that there exists a layer solution v of 

i-d..yv = fiv) znM, (2.6) 

that is, V is a solution of ( \2.6\} satisfying 

v' > m M and lim v{x) = ±1. 

Then, we have 

G'{1) = G'{-1) = (2.7) 

and 

G>G{1)=G{-1) m(-l,l). (2.8) 

(ii) Let f be any G^'''{M.) function with 7 G (0, 1). Assume that {vg}, 
with s = Sfc G (0, 1) and Sk '\ 1, is a sequence of layer solutions of 

i-dxxyvs = f{vs) ^nR, 

such that Vs{0) = 0. Then, there exits a function v such that 

lim Vg = v 

in the uniform convergence on every compact set o/R. Furthermore, 
the function v is the layer solution of 

-v" = fiv) m R 

with v(fd) = 0. 

Conditions ( 12. 7p and ( 12. Sp express that is G is of double-well type and 
/ is of bistable balanced type. Note that the statement G{1) = G{—1) 
is equivalent to 

Theorem 12.21 is actually a consequence of the following Hamiltonian 
equality and estimate, which are of independent interest. We have in- 
troduced the Hamiltonian associated to problem (II. 3p in subsection 1.1 
above. The following Hamiltonian identity for layer solutions in R 
states the conservation of the Hamiltonian in "time" — recall that x 
plays the role of time variable. Instead, the Hamiltonian inequality 
below is the analogue in dimension 1 of the classical Modica estimate 
for bounded solutions of semilinear equations Av — G'{v) = in R", 
which states that the kinetic energy is bounded at every point by the 
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potential energy, i.e., (l/2)|Vf p < G{v) everywhere in M", whenever 
G > in M. 

Notice that our Modica-type estimate is stated for n = 1. It is still 
an open problem for n >2. 

The theorem also provides an asymptotic result as s goes to 1, in 
which we recover the classical Hamiltonian identity. 

Theorem 2.3. (i) Let a E (—1, 1) and f any C^'"'(K) function, for 
some 7 > max(0,a). Let n = 1 and u be a layer solution of (11. 3p . 
Then, for every x E M. we have °° t'^\\/u{x,t)\'^dt < oo and the 
Hamiltonian equality 

(1 + a)/ -t''{ul{x,t)-ul{x,t)}dt = G{u{x,0))-G{l). (2.9) 
Jo 2 

Furthermore, for all y > and x eM. we have 
ry fa 

(1 + a) / -{ul{x,t)-ul{x,t)}dt<G{u{x,0))-G{l). (2.10) 
Jo 2 

(ii) Let f be any G^''^(M) function with 7 E (0, 1), n = 1 and {Ua}, 
with a = ak E (—1, 1) and J, —1 be a sequence of layer solutions of 
(11.31) (with u replaced by Ua for each a) such that Ma(0,0) = 0. Then, 
lima=a^4,_i Ua(-, 0) = U m the uniform G^ convergence on every compact 
set o/M, where v is the layer solution of ~v" = f{v) in M with ?J(0) = 0. 
In addition, for every x eM. we have 

lim (1 + a) / -e{ua)l{x, t)dt = -v'ix)^ = G{v{x)) - G{1) 

ai-l Jq 2 2 

and 

lim (1 + a) / -t''{ua)l{x,t)dt = 0. 

We emphasize once again that the previous estimate ( 12.10p is point- 
wise and nonlocal. 

The asymptotic result when a — )■ — 1 (i.e., s — )■ 1) of part (ii) in the 
previous theorem allows to recover from (12.91) the standard conservation 
of the Hamiltonian for the Laplacian. This will be presented in section 6 
below. 

In the case of radial solutions with limit at infinity, we establish the 
following result. Here the dimension n is arbitrary. 

Theorem 2.4. Let s E (0, 1) and f any G^'"'{M.) function, for some 
7 > max(0, 1 — 2s). 

Let n > 1 and v = v{x) = "^(1x1) be a nonconstant radial solution of 

{-Ayv = f{v) mM'^ 
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satisfying 



\x\ 



lim f (|a;|) = 0. 



oo 



Then, we have 



/(O) 



= G"(0) and G(0) > G{v{0)). 



If, in addition, v is decreasing in \x\, then 



/'(O) = -G"{0) < 0. 
The statement G{0) > G(f (0)) is equivalent to 




Jo 

As in the case of layer solutions, Theorem 12 .41 relies on the following 
statement about the Hamiltonian. In this next theorem we do not 
assume u{-,0) to have a limit at infinity. 

Theorem 2.5. Let a G (—1, 1) and f any C^'"'{M.) function, for some 
7 > max(0, a) . 

Let n > 1 and u be a bounded solution of (11. 3p which is radial in x, 
i.e., u{x,y) = u{\x\,y). 

Then, J^°° t°'\'Vu{r,t)\'^dt < oo for every r >0, and the quantity 



(1 + a) / -\ul{r,t) - ulir,t) \ dt - Giuir,0)) (2.11) 



is nonincreasing in r >0. 

3. Local realization of the fractional Laplacian and 



This section is concerned with the relation between the local problem 
(II. 3p and the nonlocal problem (II. ip . We collect also several results on 
degenerate elliptic equations with A2 weights. 

We first introduce the spaces 



where s e (0, 1) and J-' denotes Fourier transform. For Q C M""*" a 
Lipschitz domain (bounded or unbounded) and a G (—1, 1), we denote 



H\n,y'') = {u e L\n,y''dxdy) : \Vu\ e L^n^y" dx dy)} . 




RESULTS ON DEGENERATE ELLIPTIC EQUATIONS 
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3.1. Local realization of the fractional Laplacian. The fractional 
Laplacian can be defined in various ways, which we review now. It can 
be defined using Fourier transform by 

for V G iJ*(]R"). It can also be defined through the kernel representation 
(see the book by Landkof [T^ ) 

(-A)Mx) = P.V. / ]["l~jS] dx, (3.1) 

for instance for v G 5(M"), the Schwartz space of rapidly decaying 
functions. One can also define the fractional Laplacian acting on spaces 
of functions with weaker regularity. Indeed, following [18], one defines 
the space iSs(R"') of functions v such that for every k > 0, the 
quantity (1 + \x\'^~^'^'^)D'^v{x) is bounded. We denote iS^(M") its dual. 
It is then possible to check that (— A)'^ maps iS(R"') into iSs(M"). By 
duality, this allows to define the fractional Laplacian for functions in 
the space 

/:.(M"):^{.GLL(M"):£^Jg)^ 

= LL(M")n5:(R"). 

For V G CsiR"") n Cl^iW), the integral in ([31]) is well defined. This 
is clear for \x\ large. For the Cauchy principal value to be well defined 
(as X — )■ x), it suffices to assume that v is Ci^^iM.""). In particular, 
expression (13. ip defines the operator on the type of solutions that we 
consider, since they will always be bounded in M" and locally C^. See 
[iHl [H] for more comments on the various definitions of the fractional 
Laplacian and their agreement. We refer the reader to the book by 
Landkof [14] where an extensive study of integro-differential operators 
with Martin- Riesz kernels, i.e., kernels of the type (up to a normalizing 
constant) l/l^l""'"^'' is presented. 

It is well known that one can see the operator (— A)^/^ by considering 
it as the Dirichlet to Neumann operator associated to the harmonic 
extension in the half space, paying the price to add a new variable. In 
[6], Caffarelli and Silvestre proved that such a kind of realization is also 
possible for any power of the Laplacian between and 1, as follows. 

Given s G (0, 1), let a = 1 — 2s G (—1, 1). It is well known that the 
space H'{W) coincides with the trace on (9M++^ of if^(M++\ y'*). In 
particular, every v G H'^(W') is the trace of a function u G Lf^^iW^-^, y"-) 
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such that Vm G L^(R"^^, y''). In addition, the function u which mini- 
mizes 

min < / y"" {Vuf dxdy : u\Q^n+i = v \ (3.2) 
solves the Dirichlet problem 



LaU := div (y'^Vu) = in R!^+^ 
u = V on d. 



(3.3) 

n+1 ^ ' 



By standard elliptic regularity, u is smooth in M""^^. It turns out that 
—y"'Uy{-,y) converges in if^*(M"') to a distribution h G iJ^'*(R") as 
y iO. That is, u weakly solves 

div (y^Wu) = in R"+^ 

(3.4) 

-y''dyu = h on9R!^+^ 
Consider the Dirichlet to Neumann operator 

V I— J- TAv) = h := — lim y'^dyU, 

where u is the solution of (13. 3p . Then, we have: 
Theorem 3.1 ([6J). For every v G i7"(R"), 

{-AYv = dsTaiv) = -ds lim y^dyU, 

where a = 1 — 2s, dg is a positive constant depending only on s, and 
the equality holds in the distributional sense. 

In other words, given h G H^'^{W^), a function v G iJ*(R") solves 
the equation (— A)'^f = dgh in R" if and only if its extension u G 
H\Rl+\y'') solves ([32D. By duality, the same relation can be stated 
when V G £s(R"') — as it is the case of the solutions considered in this 
paper. 

3.2. Degenerate elliptic equations with A2 weights. According 
to the previous result, we must study the operator La = div (?/"V) 
in R"'*'^, whose weight y^ belongs to the class A2 since a G (—1,1). 
Through a reflection method, it will be useful to consider the equation 
in domains Q C R""*"^ not necessarily contained in R""^^. In such case, 
we extend the weight y"" by \y\°' for y < 0. That is, we define 

LaU := div {\y\^Vu) in c R"+^ (3.5) 

In a series of papers ([I0l[9]), Fabes, Jerison, Kenig, and Serapioni 
developed a systematic theory for this class of operators: existence of 



12 



XAVIER CAERE AND YANNICK SIRE 



weak solutions, Sobolev embeddings, Poincare inequality, Harnack in- 
equality, local solvability in Holder spaces, and estimates on the Green's 
function. 

In particular, as a consequence of a Poincare inequality related to 
A2 weights, they established an existence result (via the Lax-Milgram 
theorem) and a Holder continuity result. The following three results for 
La as in (13. Sp . with a E {—1,1) follow from results of [lO], stated there 
for general A2 weights. More precisely, they follow respectively from 
Theorem 2.2, Theorems 2.3.12 and 2.3.15 (and Remark 1 preceding it), 
and Lemma 2.3.5 of [lOj. 

Theorem 3.2 (Solvability in Sobolev spaces [ID]). Let n C M"+^ 
be a smooth bounded domain, h = {hi, hn+i) satisfy \h\/\y\'^ G 
L^(n, I2/I"), and g G H^{Q, \y\°'). Then, there exists a unique solution 
u G H^{^1, |?/|") of LaU = —divh in with u — g G Hq{^1, ly]""). 

Theorem 3.3 (Holder local regularity [lOj). Let C M"+^ be a smooth 
bounded domain andu a solution ofLaU = — div h in Q, where \h\/\y\°' G 
^2(n+i)^^^ |?/|"). Then, u is Holder continuous in Vt with a Holder ex- 
ponent depending only on n and a. 

Theorem 3.4 (Harnack inequality [10]). Let u be a positive solution 
of LaU = m BiR{xo) C Then, supsj,(xo)^ < CinfB^(^.Q)U 

for some constant C depending only on n and a — and in particular, 
independent of R. 

As a consequence, bounded solutions of LaU = in all of M"+^ are 
constant. 

The last statement is proved applying the previous Harnack inequal- 
ity to M — inf]Kn+i u in _B4/j(0) and letting R ^ 00. 

Corollary 3.5. Problems (13.31) and (13.41) admit at most one solution 

u with u bounded and continuous in W]^^ — up to an additive constant 
in the case of (13.41) . 

Proof. The difference w of two solutions would solve the homogeneous 
problem. We can then perform odd reflection for problem (13.31) and 
even reflection for problem (13. 4p . to obtain a bounded solution of LaW = 
in all of M"'*'^. By Theorem 13. 4[ w is constant, which finishes the 
proof. 

We remark that in case of the Neumann problem (13.41) . the previous 
Liouville and uniqueness results also follow from the Harnack inequality 
in R""^^ that we prove in Lemma [4.91 below. □ 

The existence of a bounded solution for (13.30 and (13.40 is stated 
below in this same section. 
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3.3. A duality principle. An important property of the operator La 
is the following duality property. It relates the Neumann problem for 
the operator La with the Dirichlet problem for the operator with 
the inverse weight. 

Proposition 3.6 ([6]). Assume that h G C{W), u G C'^{Mr^^), and 
y°'dyU G C(M"^"'^). If u is a classical solution of 

'div(?/" Vm) = mW\^^ 
^ = h on dWl+\ 

then w = —y°'dyU is a classical solution of 

'div(y-«Vw) = mWl^^ 
w = h on dRl-^\ 

The previous duality property is related, in dimension two, to some 
generalized Cauchy-Riemann conditions that we describe. Indeed, writ- 
ing LaU = in as 

dxiy^d^u) + dyiy^'dyu) = 0, 
we see that the associated conjugate function u is such that 



y^dyU = d^u 
-y^d^^u = dyU 



7 



hence satisfying generalized Cauchy-Riemann conditions. The function 
u is the a-conjugate of u. Similarly, u is the —a conjugate of u. Com- 
plexifying the problem by denoting = m + m, it is easy to see that 
satisfies 

dip = iy{y)dip (3.6) 
where d = + idy is the standard (9-operator and 

^ = 1 a- 

Equation (13. 6p is called conjugate Beltrami equation and has been 
extensively studied in the Calderon problem (see [2] and references 
therein) . 

3.4. Fundamental solutions. Concerning the operator La involved 
in the extension or Dirichlet problem for the fractional Laplacian, one 
has the following represention formula through a Poisson kernel. 
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Proposition 3.7 (^). Given s e (0, 1), leta=l-2s e (-1, 1). The 
function 

Ps{x, y) = pn,s- -^+27 = P. 



2;|2 _|_ ^2 J ^ ( + y"^ 



n-\-l — a 



is a solution of 

' div (y'^ VPs) = m R^^^ 
R = 6n on (9P"+^ 



(3.7) 



where 6o is the delta distribution at the origin, and pn.s is a positive 
constant depending only on n and s chosen such that, for all y > 0, 

Ps{x,y) dx = 1. 

Remark 3.8. As a consequence of Corollary 13.51 above, we have that 
for G (C n L°°)(]R"), the convolution in the x- variables 

ui-,y) = Psi-,y) * V 



is the unique solution of (Q in (C n L~)(M!^+i). 

As a consequence, for v a bounded Ci^^(W^) function, w is a solution 
of (11. II) if and only if 

u{-,y) = Ps{-,y) * V 

is a solution of (11.31) (with / replaced by (1 + a)d~^f = 2(1 — s)dj^f) 
whose trace on c^M^"^^ is v. Recall that ds is the constant from (11.41) . 

From the previous duality principle (Proposition 13. 6p and the knowl- 
edge of the fundamental solution (the Poisson kernel) of the Dirichlet 
problem (13.31) . we can find the fundamental solution of the fractional 
Laplacian or, equivalently, the fundamental solution of the Neumann 
problem (13.41) . 

Proposition 3.9 ([6j). Given s e (0,1), let a = 1 -2s e (-1,1). The 
function 

rs{x,y) = e„,,|(x,i/)|2^-" = e„,,|(x,t/)|i-'^-" 
is a solution of 

'div (|/»Vr,) = mRl+' 

-y^dyTs = 6o on dRl^^ =W, ^ ' ^ 

where 6o is the delta distribution at the origin, and en,s is a positive 
constant depending only on n and s. 
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As a consequence, rs(x, 0) = e„,s|xp'^~" is, up to a multiplicative 
constant, the fundamental solution of (—A)* in M". 

Remark 3.10. As a consequence of Corollary 13.51 above, we have that 
for h G Cc(M") {h continuous with compact support), the convolution 
in the x- variables 

= r,(-,?/) * h 

is the unique (up to an additive constant) solution of (13. 4p in {C fl 
L~)(R^). Thus, its trace 

V = * h 

is up to a multiplicative constant, the unique (up to an additive con- 
stant) continuous and bounded solution of 

{-AYv = h in M". 

Remark 3.11. The normalizing constant Cn,s in (II. 2p is given by 

-r(-s) r(i-s) r(2-s) ^ ' 

In particular, up to positive multiplicative constants, Cn,s behaves as 
s and 1 — s for s J, and s ^ 1, respectively. 

Let us also make some comments on the constant dg in the Caffarelli- 
Silvestre extension problem — see Theorem 13.11 Its value is given by 

d, = 2--^ 



r(i-.) 

The fact that this constant does not depend on n is already proved in 
section 3.2 of |6]. Its precise value appears in several papers; see e.g. 
[la [n]. Using that sr(s) = r{s + 1) and (1 - s)r(l - s) = 1(2 - s), 
we deduce, respectively, that 

ds ds 

— )■ 1 as s I and — r — )■ 1 as s t 1- (3.9) 

When u solves problem (II. 3p . its boundary condition {1 + a)du^u = 
2(1 — s)dy^u = f{u) gives that the trace of u solves 

(-A)^H,.,0)} = -A_^/(«(,-,0)). 
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4. Preliminary results on elliptic problems involving La. 

SCHAUDER ESTIMATES, MAXIMUM PRINCIPLES, AND A 
LIOUVILLE THEOREM 

This section is devoted to the proof of several general results con- 
cerning problem (11. 3p . The following definition provides the concept of 
weak solution for (11.31) . More generally, we consider the problem 

div (y" Vn) = in 5+ , , ^ , 

-y'^Uy = h on T% 

where we have used the notation introduced in the beginning of sec- 
tion 2. 

Definition 4.1. Given i? > and a function h G L^(r5j), we say that 
u is a weak solution of (14. ip if 



1/"|Vm|' G L\Bj 



and 



y^Vu-V^- / h^ = Q (4.2) 



for all ^ G CH^r) such that ^ = on r+. 

Remark 4.2. The (weak) maximum principle holds for weak solutions 
of (14.11) . More generally, if u weakly solves 

-div (y"Vn) > in 

-y^Uy > on (4.3) 

■u > on r^, 

then ti > in B^. This is proved simply multiplying the weak formu- 
lation by the negative part u~ of u. 

In addition, one has the strong maximum principle: either m = 
or u > in B^ U F^j. That u cannot vanish at an interior point 
follows from the classical strong maximum principle for strictly elliptic 
operators. That u cannot vanish at a point in Fjj follows from the 
Hopf principle that we establish below (see Proposition 14. lip or by the 
strong maximum principle of [TU] . 

Note that the same weak and strong maximum principles (and proofs) 
hold in other bounded domains of W]^^ different than 5]^. It also holds 
for the Dirichlet problem in 5^, i.e., replacing the Neumann condition 
in (USD by n > on T% 
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4.1. Schauder estimates. In the following, we prove several esti- 
mates for solutions of f 1 1.1 1) and (II. 3p . as well as for solutions of the 
Neumann problem (14. ip in B^. 

Remark 4.3. Note that the function u{x,y) = y^^"" solves LaU = 0. 
Therefore, one cannot expect the Holder regularity in y to be higher 
than than C'™'^(i'i-'=') up to the boundary {y = 0}. Thus, for s < 1/2 
(i.e. a > 0), there are solutions of LaU = vanishing on {y = 0} which 
are not up to {y = 0}. 

The next lemmas provide several regularity results. We start with 
estimates on the nonlocal equation. 

Lemma 4.4. Let f be a C^''^(R) function with 7 > max(0, 1 — 2s). 
Then, any bounded solution of 

{-Ayv = f{v) mR" 

is C^''^(R") for some < /3 < 1 depending only on s and 7. 

Furthermore, given sq > 1/2 there exists < /3 < 1 depending only 
on n. So, and 7 — and hence independent of s — such that for every 

s > So, 

Il'^llc2./3(]R") < C 

for some constant C depending only on n, so, WfWc^n, and \\v\\l^(k^) 
— and hence independent of s E (sq, 1). 

In addition, the function defined by u = Pg * v (where Pg is the 
Poisson kernel in Proposition 3.7) satisfies for every s > So, 

for some constant C independent of s E (sq, 1), indeed depending only 
on the same quantities as the previous one. 

Proof. Since v is bounded, f{y) is also bounded. Applying Proposition 
2.9 in [T8], we have 

• If s < 1/2, then for any a < 2s, w G CO'°(M"). 

• If s > 1/2, then for any a < 2s - 1, G C^^''{W). 

This implies in particular that f{v) is C"(M"). Applying now Propo- 
sition 2.8 in [T8], we have 

• If a + 2s < 1, then v G C^'^'+^'iW). 

• If a + 2s > 1, then v G Ci'"+2^-i(M"). 

Therefore, iterating the procedure a finite number of times, one gets 
that V G C^'°"(M"') for some a G (0, 1) depending only on s. Indeed, if 
q; + 2s > 1, then one can take a = a + 2s — 1. On the other hand if 
a + 2s < 1, we have that f{u) is (7°'"+^''. As a consequence, one gets 
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that u is (7°'"+*^''. Hence iterating a finite number of times, we will end 
up witli a + k2s > 1 for some integer k. This gives the C^'"" regularity. 
We now differentiate the equation to obtain 



for i = l,...,n, with fj.. and f'{v) belonging to C'^'°"(M"') provided we 
take 0" < 7. Therefore, applying Proposition 2.8 of [IH] we obtain that 
G C*'''^^^'^(M"). We iterate this procedure a finite number of times 
(as long as the Holder exponent is smaller than 7). Now, since by 
assumption 7 + 2s > 1, we finally arrive at v^^ e Ci'/3(R"), and thus 
V G C^'^(M") for some /3 > depending only on s and 7. 

For the second point of the lemma, we write the nonlocal equation 

as 



A careful look at the proof of Proposition 2.5 in [18] shows the following. 
Given sq > 1/2, if < 2 — 2so < a < 1, the operator (— A)^^'* maps 
C"(M") into C"~^"'"^''(]R") continuously with a constant independent of 
s. Here we use that the constant Cn,s in (11. 2p is uniformly bounded as 
s t 1 — see Remark 13.111 As a consequence, applying C^'^ estimates for 
Poisson equation, we deduce that v e c2'"-2+2^o(^]^n^ ^ ^^'^(M'') and 
a C^'^ estimate with a constant independent of s — indeed depending 
on the quantities in the statement of the lemma. 

We now come to the last point of the lemma. Let v G (L°° nC'^)(M"') 
for some /3 G (0, min(l, 2so)] — here we allow sq G (0, 1). We claim that 
there exists a constant C depending on n and sq, independent of s > sq, 
such that the function u defined by 



nv)v,^ in W 



n 



is (L°° n C/3)(M'^+i) with the estimate 



< C\\v\\cfi(Rn). 



Applying this fact to v, v^^ and v^^xj, we conclude the statement of 
the lemma. To prove the claim, the Poisson kernel Ps writes 




where 
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The constant pn,s is such that Hs{^) = 1 and therefore it is 
bounded uniformly in s, for s > sq. We have 



u{x,y)=l v{x - y^)Hs{Od^ 

and then 

\u{xi,yi) - m(x2, 1/2)1 



<c{\xi-X2\^ + \yi-y2f [ 



Thus we deduce the desired result taking /3 < 2so < 2s. □ 

The next lemma provides estimates for solutions of the Neumann 
problem in a half-ball. 

Lemma 4.5. Let a E (—1, 1) and R > 0. Let ip G C"{T1j^ for some 
a E (0, 1) and u E L^^^B^ji) H be a weak solution of 

-- zn B^^ C 
= 09 on Fop. 

Then, there exists /3 E (0, 1) depending only on n, a, and a, such 

that u E C^'^(bI) and y^Uy E C^'^(b^). 

Furthermore, there exist constants C\ and depending only on n, 
a, R, and also on ||^||L^(ro^) (for C}^) and ||^||c-(ro^) (for 

C^) , such that 




and 



\y''^y\\co,p(B+) - ^R- 



Proof. Multiply ip hj a cut-off function identically 1 in Fg^^g ^md call 
it if. Thus Tp is C^(M") and ^ = in Let u be the solution of 

= in Wl+^ 
= Tp on dWY^^ 

given in Remark 13.101 We have that u is continuous and bounded. 
Then, the function Mlgjj^+i solves 

{-Ayu = ds(p inW 
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and by Proposition 2.9 in [IB], we have that u G C°'^(]R") for some 
(3 G (0, 1) depending only on n and s = 

We consider now w = u — u. Then, w is a bounded function and a 
weak solution of 

'l„w = in C 

dw „ 
- — = on r!^p,„. 

Reflecting evenly the function w with respect to {y = 0}, the reflected 
function w satisfies in the weak sense the problem 

div (lyrViI;) = in 53/J/2 C M"+^ 

Since the weight |?/|" is and the function w is bounded, the regularity 
theory in (see Theorem 13.31 above) ensures that w is C^'^ for some 
/3 G (0, 1) depending only on n and s. 

Putting these two results together ensures that u G C^'^{B^ for 
some /3 G (0, 1) depending on n and s. Furthermore, we get the esti- 
mate 

for some constant C\ as in the statement of the lemma (depending on 
||v9||loo and not on ||(y9||c'<T). 

On the other hand, by Proposition 13.61 ipi^x^y) := —y°'Uy{x,y) sat- 
isfies 

L_a4j = in B^^ 

^ip = (p on r°^. 

We introduce 

i^{x,y) = Ps{-,y) * ^, 
where s is such that l — 2s = —a. Recall that Tp G C"^(M"') has compact 
support. Thus, ip is bounded and C^iMJ^-^) if /3 < min(cr, 2s) (for this, 
recall the argument on convolutions at the end of the proof of Lemma 

I43D. 

Now, we have that the odd reflection of the function ijj = ip — ip 
satisfies in weak sense 

dw{\yrV^) = in 53^/2 CM"+\ 

Hence by the results in (see Theorem 13.31 above), ip is for some 
(3 G (0,1), depending only on n and s (perhaps different than the 
previous /3). This and the above fact on ip give the desired result and 
estimates for ip = —y°'Uy = ip + ip □ 
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4.2. Gradient estimates and integrability at infinity. The fol- 
lowing two results concern bounds for solutions of problem f ll.3p . 

Proposition 4.6. Let f be a C^''^(]R) function with 7 > max(0, 1 — 2s) 
and u E L°°(R""*"^) a weak solution of problem (II. 3p . 

Then, and y"'dyU belong to L°°(M"^^). In addition, given sq > 
1/2, there exists a constant Ci depending only on n, sq, II/Hcl^ one? 
II'"IIl°o(r;;j;+1); that for every s > Sq, we have 

I|VxM|Loo(k^+1) + (1 + a)||?/"aj,M||^oo(Ign + l) < Ci. (4.4) 

Furthermore, we have 

\Vuix,y)\<— fory>0, (4.5) 

y 

where the constant C2 is uniformly bounded for a G (—1, 1). 
As a consequence of (14. 4p . we have 



Proof. The bound || VsM||j;^oo(]Kn+i) in (14. 4 p follows from Lemma and 
the fact that 

— = J- 1 as s — 1: 

2(l-s) 1 + a 

see Remark 3.11. The bound on (1 + a)\\y°-dyu\\j^oo(j^ri.+i^ in (14. 4 p follows 
from duality (see Proposition 3.6) and the boundary condition 

-(1 + a)y''dyu\y=o = f{u) G L^{W). 

The last bound ( 14. 5p follows from rescaling the equation LaU = in 
Byg/2{xo,yo) to the same equation for u{x',y') = u{xo + yox',yoy') in 
i?i/2(0, 1). Then we use that the operator La is uniformly elliptic and 
has Lipschitz coeffcients {y')"" with constants independent of a G (—1, 1) 
—since 1/2 < < 3/2 in this ball. □ 

The following result is concerned with solutions of (II. 3p with limits 
in one Euclidean variable, or in the radial variable, at infinity. 

Remark 4.7. If -u is a layer solution of (11.31) . then not only Ux{x, 0) > 
but Ux{x,y) > for every y >0. Indeed, since LaUx = in M""*"^ and 

Ux is bounded and continuous in M"^^ (by Lemma l4.4p . Remark 13.81 
gives that Ux{-,y) is the convolution oi Ps{-,y) with m^(-, 0) > 0. Hence 
the result follows. 

Lemma 4.8. (i) Let u be a bounded solution of (II. 3p such that 

lim u{x,0) = L^ (4.6) 

xi— f ±00 
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for every {x2, ■■■,Xn) G M" ^ and some constants . Then, 

f{L^) = f{L-) = (4.7) 

and 

lim u{x,y) = (4.8) 

for every {x2,---,Xn) G M"-"^ and y > 0. Moreover, for every fixed 
R> and (x2, x„) G M"""*^, ti;e have 

\\u - L^\\Lo^(B+ix,o)) ^0 asxi-^ ±oo, (4.9) 
I|Vxm|Loo(b+(^,o)) ^ as xi ^ ±cx), (4.10) 

and 

II?/XIIl-{b+(x,o)) ^ as xi ^ ±oo. (4.11) 
(ii) Let u be a radial solution of ( II. 3p s'uc/i t/iai 

lim u(|x|,0) = 0. (4.12) 

Then, 

/(O) = 0. (4.13) 
Moreover, for every fixed R> 0, we have 

\\'^\\l°°(b+{x,o)) + \\'^x'^\\l°^{b+(x,o)) + \\y'''^y\\L°°{B+{x,o)) — ^ as |x| — oo. 

(4.14) 

Proof. As in Lemma 2.4 in [S] for the half-Laplacian, the lemma follows 
easily by a compactness argument and the invariance of the problem 
under translations in Xi. Indeed, in both cases (i) and (ii) of the lemma, 
one considers the family of translated (or slided) solutions in the Xi- 
variable. In the radial case (ii) we proceed like this in each Euclidean 
variable, not only the xi. 

By the Holder estimates of Lemma H75| the translated solutions con- 
verge locally uniformly and up to subsequences, to a solution of the 
same problem f ll.3p . By assumption f l4.6p or fl4.12p . such limit is iden- 
tically constant. From this, (gZD, flA]) . fO]) . and (KW) follow 
immediately. Finally, the estimate for VxU of Lemma 14.41 leads to 

KWf or dtnD. □ 
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4.3. A Harnack inequality. The following Harnack inequality for 
linear Neumann problems will be useful in the study of stable solutions 

of (inD. 

Lemma 4.9. Let ip G he a nonnegative weak solution of 

Laifi = m B^ji C M++^ 

^ + d{x)^ = on ro^, 

where d is a bounded function in F^^. Then, 

supip < Cr inf if, (4-15) 

for some constant Cr depending only on n, a, and i?^^"||(i||£^oo(ro^). 

Proof. By scaling, one can assume R = 1. We introduce the new 
function 

The function \E'^ satisfies 

div {y^Vie'^y^'^m^)) = in S+ 

^ = ^^A{1 - a) + d{x))^^ onVl 
Therefore, choosing 

. _ IMIlL°°(rO) 

~ — 1 ' 

1 — a 

we have d^a"^^ < on F^. We consider the even extension of across 
F4, defined by 

^^(x, y) = ^^(x, -y) for (x, y) e Bi,y < 0. 

Since dr^a^"^ < on F^, satisfies 

-div (||/|"V(e-^l^l'""^^)) < in 54 

in the weak sense. 

Next, taking —A we obtain dya'^~'^ > on F^, and arguing as before 
we deduce that ^'"^ satisfies 

-div (|?/|"V(e^l^l'""§-^)) > in 54 

in the weak sense. 

Denote by Ca the operator 

Caw := div (|y|'^V(e-^l^l'"'"w)). 
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We introduce now the solutions h^^ of 

C±Ah^^ = in ^4 

These solutions are obtained from the solutions of the Dirichlet problem 
for La given by Theorem 3.2, after multiplying h"^^ by e^"^'^' °. By 
the weak maximum principle and the previous considerations, we have 
that 

< and h'^ < ^"^ in B^. (4.16) 

On the other hand, since vf^/v]/-^ = e^^*^'"" < e^^"^^"" < e^^^ in 
B^, we have that < e^^A^-A j^^^^^ ^^^^^ ^^j^a = q = 

C_Ah-^ = £A(e2^l^l'"''/i-^) and on the boundary dB^, /i^ = §^ < 
g32A^-A < g32Ag2A|2,|i-^-A ^ ^32 A A\y\^-^ f^-A^ ^^le Weak maxi- 
mum principle for the operator Ca leads to 

< e''^ je^^l^l'^/i-^} < e^^^/i-^ in ^4. (4.17) 

Next, note that La(e~"^'^'^ "/i"^) = in i?4. According to the Harnack 
inequality of Fabes-Kenig-Serapioni, Lemma 2.3.5 of [10] (Theorem 3.4 
above), we deduce that 

sup(e-^l^l"''/i^) < Cinf(e-^l^l"''/i^) 

Bi 

for some constant C depending only on n and a. Thus, 

sup/i^ < Ce^inf /i^. (4.18) 
B+ B+ 

Using (??) and fl4.17p . we deduce 

< < /i^ < e^^^/i-^ < e^^^§-^ < e'^^^<^ in 5+ (4.19) 

Finally, f l4.19p and f l4.18p lead immediately to the desired result. □ 

4.4. A Liouville theorem. We prove the following theorem, which 
will be useful in [4J to prove a symmetry result for stable solutions of 
(11. 3p in . This is a generalization to degenerate elliptic equations 
of the Liouville theorem given in fSj. This type of result had been 
already used in [1] to prove the De Giorgi conjecture for reactions in 
the interior in three dimensional spaces. 
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Theorem 4.10. Let ip G Lj^j,(M" ) he a positive function. Suppose 
that a e Hl^^{W^^^ .y"") is such that 



-(Tdiv (y V^Vct) < m 
-ay^dya < on dRl^^ 



(4.20) 



in the weak sense. Assume that for every R> 1, 

y\avf dxdy < CR^ (4.21) 



for some constant C independent of R. 
Then a is constant. 

Proof. We adapt the proof given in [Sj. Let C be a function on 
[0, +00) such that < C < 1 and 



c 



1 for < t < 1 
for t > 2. 



For R> I and {x,y) E R++\ let CR{x,y) = C{r/R), where r = \{x,y)\. 
Muhiplying f l4.20p by ("^ and integrating by parts in M""*"^, we obtain 



yXl <^'l Va|2 dxdy < -2 / y^R ^'"cr VCr Va dxdy 



jn + l 



< 2 



n + l 



yXW\Va\^ dxdy 



n + l 



r\{R<r<2R} 



1/2 



yVV^lVCfll" dxdy 



n + l 



1/2 



< c 



yXl^^\Va\^ dxdy 



n + l 



r\{R<r<2R} 



1/2 



— / ?/"(^o-) dxdy 



1/2 
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for some constant C independent of R. Using hypothesis (I4.2ip . we 
infer that 



/ 



yXl^'\Va\' dxdy< 



n+l 



c 



1 1/2 



yXW\ya\' dxdy 



^r\{R<r<2R} 



(4.22) 



again with C independent of R. Hence, /j^n+i V^^l Vap dxdy < C 

and, letting i? — )■ oo, we deduce /^n+i y"''{P'\V a\'^ dxdy < C. It follows 

that the right hand side of fl4.22p tends to zero as i? — ?■ oo, and therefore 
Jj^„+i Vcrp dxdy = 0. We conclude that cr is constant. □ 



4.5. A Hopf principle. The following proposition provides a Hopf 
boundary lemma in our context. 

Proposition 4.11. Let a G (—1, 1) and consider the cylinder Cr^i = 
r^j X (0, 1) C M""^^ where T^j is the ball of center and radius R inW. 
Let u e C(C^) n H\CR^^,y'') satisfy 

LaU < in Cr^i 
M > in Cr^i 

m(0,0) = 0. 

Then, 

„u(0,y) 
limsup -y" ^ < 0. 

y^o+ y 

In addition, if y"-Uy G C^Cr^i), then 

d^au{0,0) < 0. 

Proof. Consider the function on Cr^i defined by 

WAix, y) = y^^iy + Ay'^)ip{x), 

where A is a constant to be chosen later and = ^{x) is the first 
eigenfunction of — in r^y2 with Dirichlet boundary conditions, i.e., 

-A^.(/3 = \i(f in r°j/2 
^ = on dT%^. 
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Notice that Ai > and that we can choose ip > in r^y2 with — 
1. The function wa satisfies 

LaWA = ^{x) 1^(2 -a)-Xi{y + Ay^) | in Cr/2,i 



Wa > in Cr/2,1 

wa = o on X [0, 1). 

Therefore, choosing A large enough, we have in Cr/2,i 

LaWA > 0. 

Hence, for £ > 0, 

La{u - EWa) < in Cr/2,1 

and u — EWa — u > on dT^^^^ x [0, 1). Moreover, taking £ > small 
enough, we have on F^^g ^ {v — 1/2} 

u > EWa, 

since u is continuous and positive on the closure of this set. Notice 
furthermore that wa — on T^j^ x {y — 0}. Thus, we have 

La{u - SWa) < in Cr/2,1/2 

u - ewA > on dCR/2,i/2- 



The weak maximum principle then implies that in Cr/2,i/2 

u — SWa > 0. 

Consequently, this leads to 

limsup-y«^^^^ < £limsup-y"^^^^^^ = -£(^(0) < 0, 

y^o+ y 2/-J-0+ y 

as claimed in the proposition. 

Assume, in addition, y^Uy G C{Cr^i). Let yo < 1/2. Since {u — 
£Wa)(0, •) > in [0, yo] and {u — ewA){0,0) — 0, we have {uy — 
e{wA)y){0,yi) > for some yi e (0, yo)- Repeating this argument for 
a sequence of ?/qS tending to 0, we conclude that ~y°'Uy < —ey"'{wA)y 
at a sequence of points (0,?/^) with yj J, 0. Since we assume y"'Uy con- 
tinuous up to {y = 0} and —e{y"'{wA)y){0,yj) —eip{0), we conclude 
that d^au{0, 0) < 0. □ 

Corollary 4.12. Let a G (—1, 1) and £ > 0. Let d be a Holder contin- 
uous function in r° and u G L°°{B^) fl H^{B^, y") be a weak solution 



28 XAVIER CAERE AND YANNICK SIRE 

of 



LaU = in 
M > in B'^ 

d„au + d{x)u = on 
Then, u > in B^ U r° unless u = in Bf. 

Proof. We apply Lemma 14.51 to obtain that u (for this, see the proof of 
the lemma) and y°'Uy are up to the boundary. Hence the equation 

d^au + d{x)u = (4.23) 

is satisfied pointwise on r°. If m is not identically in Bf then u > 
in Bf by the strong maximum principle for the operator La- Now, 
if u(a;o,0) = at some point (xo,0) G r°, then a rescaled version of 
Proposition 14.111 gives dyau{xQ,Q) < 0. This contradicts fl4.23p . □ 

4.6. A maximum principle. Here we present a maximum principle 
related to the operator La and to the fractional Laplacian. We will 
need it in our subsequent article to prove monotonicity properties for 
solutions in M with limits, as well as the uniqueness (up to translations) 
of layer solutions in R. Recall that section 3 already contained some 
Liouville and maximum principles for these operators. 



Lemma 4.13. Let m G (C n L°°)(R^+i) with y^Uy G C(M!^+i) satisfy 

LaU = m Wl^^ 

dv , (4.24) 

^ + d(x)u>0 on dRl+\ ^ ' 

where d is a hounded function, and also 

m(x,0)->0 as \x\^oo. (4.25) 

Assume that there exists a nonempty set H C such that u{x, 0) > 
for X & H , and d{x) > for x ^ H. 

Then, u > %n W^. 

Proof. By Remark [3.81 applied to v{x) = u{x, 0) — infjjn «(-, 0) > 0, we 
see that u — infjRn u{-, 0) > in W^-^. Thus, infj^n+i u = infRn u{-, 0). 

Arguing by contradiction, assume that there exists a point (xq, yo) in 
M"^"*^ such that u{xo, yo) < 0. Then, in case infj^n+i u = 0, the minimum 
of M is a achieved at (a;o, yo)- In case infjjn+i u = inf^n «(-, 0) < 0, using 
that u{x,0) — )■ as |x| — )■ +oo, there exists a point (xi,0) at which 
the minimum of u is achieved. In both cases we conclude that the 
nonpositive minimum of u is achieved at a point {x2,y2)- 
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By the strong maximum principle, we cannot have 1/2 > 0, since u 
is not identically constant (recall u(-,0) > in if 7^ 0). Thus 1/2 = 0. 

According to the Hopf lemma 14.111 and since y"'Uy G C(M';+^), we have 

du 



(0,Z/2) <0. 



Since u(a;2,0) < then X2 4- ^1 thus we have ^(0:2) > 0. Now, 
using the boundary condition in f l4.24p at x = X2, we reach a contra- 
diction. □ 

Remark 4.14. Lemma [4.131 can be stated in an equivalent way using 
the equation 

(-A)'t; + d{x)v > in 

and assuming the same conditions on v as those for 0) in the pre- 
vious lemma. In addition, an alternative proof of the lemma can be 
given using the integral expression (11. 2p for (— A)*f (a;2), that will be 
negative at a point of minimum (since v is not identically constant in 
the proof). 



5. HAMILTONIAN ESTIMATES 

This section is devoted to establish the main facts needed to prove 
Theorems 12. 3l and l2. 51 We start with an easy lemma that will be needed 
later in several occasions. 

Lemma 5.1. Let u G L°^{W^^^ he a bounded solution of (11.31) . Then, 
for all X G R", we have f^°° t"'\Vu{x,t)\'^dt < 00. In addition, the 
integral can be differentiated with respect to x E R" under the integral 
sign. Furthermore, 

r+oo 

lim / t''\Vu{x,t)\'^dt = (5.1) 



M 



uniformly in x E R". 

// in addition, u is either a layer solution in R (here n = 1) or u is 
a radial solution in M" for which lim|^|_^oo w(|x|, 0) exists, then 

r+oo 

lim / r|VM(x,t)|2rft = 0. (5.2) 

I^^Koo Jo 

Proof. The first two statements and (15. ip follow directly from the gra- 
dient bounds in Proposition 14. 6[ The statement (15. 2p is a consequence 
of (ISTT!) and of Lemma Sill □ 
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5.1. Hamiltonian equality and estimate for layer solutions. This 
subsection contains two lemmas. Ttie first one establislies tliat the 
Hamiltonian is conserved for layer solutions in dimension one. 

Lemma 5.2. Let n = 1 and assume that u is a layer solution of (11. 3p . 
Then, for a// x G M we have °° t"'\Vu{x,t)\'^dt < oo and the Hamil- 
tonian identity 

p+oo j.a 

(1 + a) -{ul{x,t)-ul{x,t)}dt = G{u{x,0))-G{l). (5.3) 
As a consequence, 

G(1)=G(-1). (5.4) 

Proof. The integrability of t°'\'Vu{x,t)\'^ follows from Lemma 5.1. We 
now establish equality (15.31) . It will be crucial that the weight in La 
does not depend on the tangential variable x. 
Following [S], we consider the function 

— {ul{x,t) - ul{x,t)} dt. (5.5) 

Lemma 15.11 allows us to differentiate under the integral in (15.51) to get 
d 



^^V{x) = J ^{UxxUx - Uj:yUy){X,t)dt. 

Noticing that 

LaU = dy{y'^Uy) + y^u^^x = 
and after an integration by parts (which is justified by Lemma [5. ip we 
have 

-^v{x) = lim y''uy{x,y)ux{x,y) = --^—-^G{u{x,0)). 

dx y~*o+ I + adx 

The function (1 + a)t>(x) — {G{u{x, 0)) —G{1)} is then constant in x. 
Letting x — )■ +oo and using Lemma 15.11 we have that this constant is 
actually zero. Letting now x — > — oo and using Lemma ISTTl we deduce 
G(1)=G(-1). □ 

We have obtained that a necessary condition for the existence of a 
layer solution in M is that G{1) = G{—1). The other necessary condi- 
tion will follow from the following result — our Modica-type estimate 
for layer solutions in dimension 1 (Theorem 12. 3p . 

Lemma 5.3. Let n = 1 and assume that u is a layer solution of (II. 3p . 
Then, for every y > and all x we have 

(1 + a) r^-{ul{x,t)-ul{x,t)} dt<G{u{x,0))-G{l). (5.6) 
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Proof. We introduce the function 

w{x,y) = J — {ul{x,t) ~ul{x,t)] dt, 

which is bounded in all by Lemma 15.11 We introduce the function 

w{x,y) = -^\g{u{x,0)) -G{1)] -w{x,y). 

The function W is bounded in and we need to show that w > in 

We first derive some equations for W which will be useful in the 
sequel. We have, for all y > 0, 

Wy{x, y) = -y (M^(a;, y) - u\{x, y)). (5.7) 

Furthermore, using LaU = and integrating by parts as in the previous 
proof, one gets for all ?/ > 

y) = y^u^ix, y)uy{x, y). (5.8) 

Using the two previous equalities and the equation LaU = 0, we have 
for all ?/ > 

LaW = -ay'^'^-^l (5.9) 

and 

L_aW = -ay'^u^y. (5.10) 

We claim that w does not achieve its infimum at a point in M^. We 
assume the contrary and reach a contradiction. Let (xo,i/o) be a point 
where the infimum is achieved. There are now two cases depending if 
I/O is on the boundary or not. We will also use that w is not identically 
constant. Indeed, if it were, since w(-, 0) = then 

constant = W(-,0) = ^^|g'(u(-, 0)) - G(l)}. 

Thus G is constant in (—1, 1), / = in (—1, 1) and m is a bounded 
function satisfying (11. 3p with / = 0. Hence, after an even reflection 
across {y = 0}, Theorem 3.4 ensures that m is a constant, a contradition 
with > 0. 

Case 1: yo = 0. After a translation in x, we may assume a;o = 0. 
Since = is a global minimum of (1 + a)w{-, 0) = G{u{-, 0)) — G{1), 
we have 

= {d/dx)G{uix, 0))U=o = -/(m(0, 0))m,(0, 0), 
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and therefore 

= -/(M(0,0)) = (l + a) lim yX(0,l/), (5.11) 

since m is a layer solution (i.e., Ux{x, 0) > 0). For every (x, y) G M^, we 
have by Remark 14. 7[ 

Ux{x,y) > in M^. (5.12) 

We now divide the conclusion into two subcases. Let consider first 
the case a > 0. By (15. 7p . we see that y"'Wy is Holder continuous up 
to y = 0. Since LaW < by (15. 9p and w is not identically a constant, 
we have w > w{0, 0) in M^. Thus the Hopf principle (see Proposition 
14. lip gives that 

> - lim y^WyiO^y). 

Now, using (15. 7p and (15. lip , we have 

> - lim y^Wy{0,y) 

y2a 

= lim—{ul{0,y)-ul{0,y)} 
j/-i-o+ z 

2a 

= \imJ-ul{Q,y)>Q, 

a contradiction. 

We turn now to the case a < 0. Since (0, 0) is a global minimum for 

W{x, y), one gets 

> liminf -y"''Wy(0,?/) 

J/-J-0+ 

= lim inf ^ {uliO, y) - 7/^(0, y)) = ]^ul{Q, 0) > 0, 

a contradiction. We have used that, by Lemma 1475) [^.^(O. y)\ < Cy^"" — )■ 
as ?/ 0+. 

Case 2: yo > 0. By (I5.12p . we have > in M^. Using (15. 8 p and 
(15.1 op . we obtain 

= + ay~'^ul = L_aW + {ay~^~"'—)wx 

which is the same as 

= V ■ (y'^Vw) + b{x, y)wx inM^, 

with b{x, y) := ay~^~°-UyU~^ . But this last operator is uniformly elliptic 
with continuous coefficients in compact sets of > 0}. Thus it cannot 
achieve its minimum at (xo,?/o), since ?/o > and we have proved that 
w is not identically constant. 
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Therefore, we now know that w cannot achieve its infimum at a point 
in R^. To finish the proof, assume first 

inf w < 0. 

By Lemma 15.21 w{x, y) — > as ?/ — )■ +oo locally uniformly in x. By 
Lemma we have w{x,y) — )■ as |x| — > +oo uniformly in y. There- 
fore, the infimum of w being negative, it should be achieved at a point 
in R^, a contradiction with what we have proven. Therefore, 

inf w > 0, 

1II>2 



Thus, if w vanished at some point in Rl, this point would achieve 



i.e. w > 0. _ 

the infimum of w, a contradiction. Hence w > in R^ as stated in the 
lemma. □ 

5.2. The Hamiltonian for radial solutions. The next lemma deals 
with bounded radial solutions u of (11.31) . Here we do not assume u to 
have a limit at infinity. 

Lemma 5.4. Let u be a bounded solution of i\1.3\\ . Assume that u = 
u{\x\,y) is radially symmetric in x. Then, 

hoo 

is a nonincreasing function ofr. 
Proof. The function u solves the 

Urr + + Uyy + ^Uy = in (0, +Oo) X (0, +Oo) 

-(1 + a)?/% = /(m) on (0,+cx)) X = 0}. 

Let 

-[ul{r,t)-ul{r,t)\dt. 



(1 + a) / ^-\ui{r,t)-ui.{r,t)\dt-G{u{r,Q)) (5.13) 







By Lemma [5.11 we can differentiate under the integral with respect to 
r and obtain 

dw{r) '•+°° 
dr 

Performing one integration by parts and using the equation and Lemma 
15. H we end up with 



dw{r) _ n 1 /""^"^a^2(^^^^^^^^LG"(M(r,0))M,(r,0). 



dr 1^ .10 ^ 1 + a 
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As a consequence, the function 

{l + a)w{r) ~G{u{r,0)) 
is nonincreasing in r, as claimed. Furthermore, 



oo 



^{{l + a)w{r)-G{u{r,0))} = -{l + a)'^ I t-ui{r,t) dt. 
dr r 



(5.14) 

□ 



6. The limit s 1 and the classical Laplacian 



In the following, we investigate the asymptotic s — )■ 1. For this, we 
will use crucially the previous Hamiltonian estimates. We prove the 
following theorem. 

Theorem 6.1. Assume that f G for some 7 G (0, 1) and that 

{vs}, with s = Sfc G (0, 1) and Sk '\ 1, is a sequence of layer solutions 
of 

{-dxxYvs = f{vs) in M, 
such that fs(0) = 0. Then, there exits a function v such that 

lim Vs = v 

in the uniform convergence on every compact set ofW. 
Furthermore, the function v is the layer solution of 

-v" = f{v) m M 

with v{Q) = 0, and satisfies the Hamiltonian equality 

^(v'Y = G(v) - G{1) inR. (6.1) 

The previous theorem is stronger than just saying that the limit 
when s goes to 1 is a solution of an ODE, since it states that the limit 
is actually a layer itself. We can see Theorem 16.11 as a stability result 
in the class of layer solutions of nonlocal (and local) equations. 

Proof of Theorem \6.1[ Let Vg be a layer solution of 

{-d^j^xYvs = fivs) in R 

with Vs{0) = 0. Then, the extension Ua = Ps * Vg of Vs satisfies 

div {y^Vua) = in 

I dUa 



:i + a)^ = Caf{ua) onM^ 



on dRl, 
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where Ca = dj^{l + a) = d~^2{l — s) and ds is the constant in Theorem 
13.11 and Remark 13.111 By f l3.9p , we know that tends to 1 as a goes 
to —1. The weak formulation of this problem is 

(l + a) [ y'^Vua ■ Ve - / CafiUa)^ = (6.2) 
Jr\ Jr 

for all ^ G C^(]R+) compactly supported. 

First notice that, by the regularity result in Lemma which is uni- 
form as s t 1, the functions Ua and d^Ua converge over compact sets (up 
to a subsequence) to a function u_i = u^i{x,y) and its x— derivative 
as a — )■ — 1 (which corresponds to s — )■ 1). We now choose the following 
test function: ^{x,y) = rii{x)ri2{y) , where ri2{y) = 1 for < ?/ < 1 and 
V2{y) = for ?/ > 2, whereas rji is any test function. We deduce 



:i + a) 



y'^^r][{x)r]2{y)dxUa + r]i{x)ri2{y)dyUa^dxdy 



Caf{ua)r]i{x)dx = 0. 

We now pass to the limit in each term. Thanks to the uniform 
bounds of Lemma 14.41 we have that 



lim Ca / f{ua)rii{x) dx 

f{u-i{x, 0))rii{x) dx. 

Note that the measure (1 + a)y"' dy is a probability measure on (0, 1) 
converging as a J, —1 (in the weak—* sense of measures) to the Dirac 
measure ^o- More precisely, given functions Wa = Wa{y) continuous in 
[0, oo), with \yWa{y)\ < C in [0,+oo) (with C uniform in a) and with 
Wa converging to a function W-i uniformly in compact sets of [0, +oo), 
then 

POO 

lim (1 + a) / y''wa{y)dy = W-i{0). 

ai-l Jq 

Indeed, given e > let 6 > such that \wa{x) — Wa{0)\ < e for all 
X G (0,(5). Then, we write 

/•oo 1*5 foo 

(1 + a) / y^Waiy)dy= (l + a) y''Wa{y)dy+{l + a) y^Waiy)dy. 



We have that 
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(1 + a) / y''iw^i{0)+WaiO)-w^i{0))dy 
Jo 

tends to lima^_i w_i(0)(5^+'^ = w_i(0) and 

(1 + a)/ y''\wa{y)-wamdy<6'+^e. 
Jo 

Finally, 

(1 + a) / y-\wa{y)\dy<{l + a)C y-Uy < C^—^S'' 
Js Js -a 

(6.3) 

as a J, —1. This proves the claim above. 
We now divide the integral 



= (1 + a)/ y"-r][{x)r]2{y)d^Uadxdy 

JMx(0,+oo) 

+ (l + a)/ y''r]i{x)r]'^{y)dyUadxdy. 

Jrx(1,2) 

Thanks once again to Lemma 14.41 the observation above (with Wa 
V2{-)dxUa{x, ■)) and the gradient bounds of Lemma HSl we deduce 

lim(l + a) / y''ri[{x)r]2{y)dxUadxdy 



x(0,+oo) 

u'_i{x, 0)ri[{x) dx. 

By the same lemma, |5j^Ma| < Cy^^ uniformly in a G (—1, 0), and thus 
the same computation as in (16. 3p shows that 

lim (1 + a) / y''T]i{x)r]2{y)dyUadxdy = 0. 

"^"^ jRx(l,2) 



lx(l,2) 

Therefore, the function v := v{x) = u^i{x, 0) satisfies 



^^71 - fiv)r]i = 0. (6.4) 
Hence U is a weak solution of 

-v" = f{v) in M, 
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such that U(0) = and v' > in M. As a consequence, the function v 
admits hmits at ±oo, 

hm v{x) = G [-1, 1]. 

We now prove the convergence of the Hamiltonian, which will provide 
in addition that the function v is actually a layer, i.e. L"^ = ±1. We 
apply the Hamiltonian estimate (15. 6p with y = to the layer Ua for 
some a G (— 1, 1) with G replaced by CaG. We deduce 

<G-G{1) in (-1,1). (6.5) 

Next, we integrate the equation satisfied by v, we use the above 
observation now with Wa{-) = {dxUa)'^{x, ■) and we use the Hamiltonian 
identity (15. 3p for the layer Ua to obtain for all x G M 

G{v{x)) - G{L^) = -(U')'(x) = lim / y'^{dxUa)\x,y)dy 

(6.6) 

= lini|il±^ / °^y''idyUa)''ix,y)dy + cMuaix,0))-CaGil)], 

ai-l I Z Jq J 

and thus 

G{v{x))-G{L+) > limCa{G{ua{x,0))-G{l)) = G{v{x))-G{l). (6.7) 

a^— 1 

Hence we have that 

G(L+) < G(l), 

that together with (16. 5 p and > (since v{0) = 0) gives = 1. 

In addition, we deduce that the inequality (16.71) must be an equality. 
Thus, the term that we have dropped to obtain the inequality must be 
zero, i.e. 

lim y\dyUaf{x,y) dy = 0. (6.8) 

a->-l 2 Jq 

In the same way, we prove that = —1. Hence v is the layer solution 
connecting —1 to 1, with v{0) = 0. The uniqueness of such v follows 
from the Hamiltonian equality (16. ip . □ 

7. Proof of Theorems EIS EH, and 

We prove in this section the main theorems of our paper. They will 
follow easily from our results in previous sections. 

Proof of Theorem \2.3[ Part (i) follows from Lemmas 15.21 and 15.31 Part 
(ii) follows from Theorem 16 . 1 1 and from (16. 6p and (16. 8 P in its proof. □ 
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Proof of Theorem \2.<X If w is a layer solution of fl2.5p . its extension u is 
a layer solution of fll.3p in M^, up to a multiplicative constant in front 
of / that tends to 1 as s ^ 1- In part (i) of the theorem, (12.71) follows 
from (14. 7p in Lemma 14. 8[ The equality in (12. 8p is (15. 4p of Lemma 15. 2[ 
while the inequality in (12. 8 p follows from taking y = in the statement 
of Lemma 15. 3[ Part (ii) of the theorem follows from Theorem 16.11 □ 

Proof of Theorem \2.5[ It follows from Lemmas 15.11 and 15. 4[ □ 

Proof of Theorem \2.4\ If f is a radial solution in M", its extension 
u = u{x,y) belongs to L°°{W^^) and it is a solution of (II. 3p . up to 
a multiplicative positive constant in front of /. Clearly, m is a radial 
solution in X G M". 

The relation /(O) = is (KWf in Lemma |Ml 

The conclusion G{0) > G{u{0, 0)) of the theorem follows from Lemma 
15.41 and (I5.14p . Indeed, in (15.131) we let r = and later r ^ oo, and we 
use (15. 2p in Lemma [5. II to obtain 

-G{v{{))) > -(1 + a) / -ul{0, t) dt - G{v{0)) > -G{0). 



2 



ThusG(O) > G(t;(0)).But ifG(O) = G(t;(0)) then the function in ([513]) 
would be constant in r. Hence, by (I5.14p and since n > 1, Ur{r,t) = 
for all r and t, and then u is constant, contrary to our assumption. 
It only remains to prove the other statement of the theorem, 

/'(O) = -G"{0) < 

under the assumption Ur < 0. Without loss of generality, and to sim- 
plify notation, we may replace / by a positive multiple of / in ( II. 3p 
and hence assume that u solves 

nil 

— = fiu) on 9R"+i = M". 

We differentiate both equations with respect to r = using that the 
first one reads Urr + '^^Ur + Uyy + ^Uy = in (0, +oo) x (0, +oo). Let 

ilj:=-Ur>0 in (M" \ {0}) X (0,+oo). 
We deduce that 

71 — 1 

div (y" VV^) = ^rl/> in (»" \ {0}) x (0, +oo) 

and 

f'(u)ih = for y = 0. 
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For Xo G M", let ip^°{x,y) := ip{x — Xo,y), a positive function in 
(M" \ {xo}) X (0, +00). Let m''°(x, y) := u{x — Xq, y). We have 

11 — 1 

div iy^ W'^o) = 1 r72/>^° in (»" \ {xo}) x (0, +00) (7.1) 

\x — Xo| 

and 

f\^^')r' = for y = 0. (7.2) 

For i? > 0, consider the cyhnder Cr = T^^^ (0, R) C M^+\ where 1°^ 
is the ball of center and radius R in R". Let ^ be any function in 
Cr vanishing on {|x| = i?} x [0,-R) and on F^j x {?/ = R}. 

For |xo| > -R, we multiply (17.21) by ^'^/ip'^° — note that ip''^ > in 
F^j X [0, R) — , we integrate in F^j. and use (17.11) to obtain 



^ R 



[ div(y«VV^^°)^+i/'^VV^^°-V^ 

Jcr i' ° ^ " 

2 



n 


- 1 


\x - 


1 9 

- a;o| 


n 


- 1 


\x - 


- Xop 


n 


- 1 



y ^ +y i^^v^ ^ ^^^^^2 



Jcr k-xor 

where we have used Cauchy-Schwarz inequality. 
Letting |xo| — )■ 00, we deduce 



/'(O) < \, (7.3) 



/roe 



^2 

R 



for any function ^ in C/j vanishing on {|x| = R} x [0, i?) and on 
Tl^{y = R]. 

Let = ^Pr{x) > be the first eigenfunction of — in F^j with 
Dirichlet boundary conditions, i.e., 



-A^Lpn = in Vj^ 



yDR = on 9F° 



R 

I 

R^ 



where Xr = c{n) /i?^ > is the first Dirichlet eigenvalue of — A^. in the 
ball F^j. Let hn = h^i^y) G [0, 1] be a smooth function with compact 
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support in [0, R) such that /i^j = 1 in [0, R/2] and \h'j^\ < C /Rioi some 
constant C. Take 

in (O). We have 

Jcr Jcr 

pR rR 

\r 



[ y''hldy+ [ y^ih'^ydy] [ ^1 
Jo Jo J Jr% 

CR-^ j\''hldy + j\\h'^fdy^ j^^ i 

rR r p2-2s r 

< CR-' / y'^dy / eR = CR-^- / 

Jo Jr% 2 -2s 



R 

R 



^ R-'' I eR. 



2 -2s 



^ R 



Using this in fl7.3p and letting i? — )■ oo, we conclude that /'(O) < 0, as 
claimed. □ 
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